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Abstract. One of the goals of COMPASS experiment [I] is the determination of the gluon 
polarisation, AG/G, for a deep understanding of the spin structure of the nucleon. The gluon 
polarisation can be measured via the Photon-Gluon-Fusion (PGF) process. One of the methods 
^ \ to identify this process is selecting high pr hadron pairs in the final state [2 . The data used for 

this analysis were collected by the COMPASS experiment during the years 2002 to 2006, using 
O • a 160 GeV naturally polarised positive muon beam impinging on a polarised nucleon target. A 

D I new result of AG/G from high pr hadron pairs in events with > 1 (GeV/c)^ is presented. 

■ This result has a better precision due to the addition of 2006 data and an improved analysis 
based on a neural network approach. The gluon polarisation is also presented in three bins of 

> 
o 

I 1. The Gluon Polarisation and The High pT Analysis Formalism 

CN ■ The nucleon spin sum rule can be written in a heuristic way as: | = + + L, where AS 

and AG are, respectively the quark and gluon contributions to the nucleon spin and L is the 
orbital angular momentum contribution coming from the partons. The purpose of this analysis 
is to estimate the gluon polarisation, AG/G. The analysis is performed in two complementary 
kinematic regions: < 1 (GeV/c)^ (low Q^) [3J and > 1 (GeV/c)^ (high Q^) regions. The 

■ present work is mainly focused on the analysis for high Q^. 

^ . Spin-dependent effects can be measured experimentally using the helicity asymmetry 

defined as where (^) and (^) refer to the parallel and anti-parallel spin helicity 

configuration of the beam lepton (^) with respect to the target nucleon (<^ or According to 
the factorisation theorem in DIS, the (polarised) cross sections can be written as the convolution 
of the (polarised) parton distribution function, (A)g^, the hard scattering partonic cross section, 
(A)(j, and the fragmentation function. 

The gluon polarisation is measured directly via the Photon-Gluon Fusion (PGF) process, 
which allows to probe the spin of the gluon inside the nucleon. To tag this process directly in 
DIS a high pT hadron pair data sample is used to calculate the helicity asymmetry. Two other 
processes compete with the PGF process in leading order QCD approximation, namely the 
virtual photo-absorption leading order process (LP) and the gluon radiation (QCD Compton) 
process. The helicity asymmetry for the high pT hadron pair data sample can thus be 
schematically written as: 
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Ai^ixB,) = RPGF a{:r -^(xg) + Rlp D A^'^ixBj) + Rqcdc ai;^L Ai^"(xc) (1) 

The Ri are the fractions of each process, i refers to the different processes. a^L represents 
the partonic cross section asymmetries, AaYa^ also known as analysing power. D is the 
depolarisation factor, which is the fraction of the muon beam polarisation transferred to the 

virtual photon. The virtual photon asymmetry is defined as = ^-^^2^- This 

asymmetry is estimated using a parametrisation based on the Ai asymmetry of the inclusive 
data [4J. A similar equation to ^ can be written to express the inclusive asymmetry of a data 
sample, Af^^K 

Using equation ^ for the high pT hadron pair sample and the above mentioned equation for 
the inclusive sample the following expression is obtained: 



A = ai - a2 , A'°'' = -A,{xBj)D^ - Ai{xc)f3i + Ai(x^)/32, 
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Due to the fact that AG/G is present in formula (Ej) at two different xq (noted xq and x^), 
the extraction of AG/G requires a definition of the averaged Xq^ at which the measurement is 
performed. 
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2. Data Selection 

Data from years 2002 to 2006 are used. These data were obtained from polarised muons of 190 
GeV/c scattered off a polarised LiD target at the COMPASS experiment at the CERN SPS. The 
selected events have an interaction vertex containing an incoming muon and a scattered muon. 
The data are divided into two sets: the high pT hadron pair and the inclusive data sample. 
Both data sets have the > 1 (GeV/c)^ kinematic cut applied. Another cut is applied on the 
fraction of energy taken by the virtual photon, y: 0.1 < y < 0.9. These cuts are used to select 
the inclusive sample. 

In the high pT hadrons data sample, besides the inclusive selection, at least two outgoing 
high pt hadrons are required. These so-called hadron candidates must fulfil the following 
requirements: the hadrons with the highest transverse momentum must have pT > 0.7 GeV/c 
and the hadron with second highest transverse momentum, pT > 0.4 GeV/c . This requirement 
constitutes the high pT cut. 



3. The Weighted Method Approach and The Neural Network for AG/G Extraction 
The purpose of this analysis is to calculate the gluon polarisation, AG/G, in an event-by-event 
basis using an optimal weight which improves the figure of merit. The asymmetry used for the 
AG/G extraction is related with the experimental helicity asymmetry, ^l^, described in section 
dl using a weigthing factor w. The correct weight, in this case, should take into account all the 



variables that appear in the set of equations ([3]) to ([5]), namely: w = fDPi,X^ where / is the 
dilution factor, the fraction of polarisable target material, D is the depolarisation factor, the 
fraction of muon polarisation is transferred to the virtual photon, P5 is the muon polarisation 
and A is defined in equation [3l 

In this analysis it is not possible to tag the events of each involved process, therefore the 
process fractions Rj and the partonic asymmetries a^L cannot be directly determined from the 
data samples. To estimate or parametrise these quantities a neural network [S] is used. The 
neural network is trained by taking as input samples obtained by Monte Carlo (MC) simulation. 
The result is a parametrisation which is used for the real data to provide the values for all these 
quantities in an event by event basis. 

4. MC Simulation 

The event simulation is one of the issues of major importance for this analysis since several 
parameters for the AG/G extraction are taken from the simulation. Thus a strong effort was 
made to achieve a simulation very close to the real data. The MC production comprises three 
steps: first the events are generated, then the particles pass through a simulated spectrometer 
using a program based on GEANT 3 [6J and finally the events are reconstructed using the same 
procedure applied to real data. For the first step the LEPTO 6.5 [7] DIS event generator is used 
together with a leading order parametrisation of the parton distributions. The fragmentation 
is based on the Lund string model [8] implemented in JETSET [9j . Figure [T] illustrates the 
quality of the simulation obtained at the end. 

The next step was the tuning of the event generator. The Parton Distribution Function 
(PDF) set used in this analysis is MSTW2008LO (H]. For the calculation of Fl in LEPTO the 
ratio R{x,Q^) = ^ is needed. Here the parametrisation of ref. [7J was used. To improve the 
description of the hadrons, higher order QCD corrections are partially simulated by including 
gluon radiation in the initial and final states {parton shower - PS) [7]. An extensive and 
careful tuning of the Lund string fragmentation function parameters and of the hadron intrinsic 
transverse momentum parameters was performed. 

The tuning was done first on the fragmentation function parameters, in this case of the the 

Lund string function: f{z) oc ^(1 — z)^exp{— ^ ^-^ ) where rri]_ = m? the parameters are a 

and b which are controlled by the JETSET parameters PARJ(41) and PARJ(42). 

The second step was the tuning of the hadron intrinsic transverse momentum model, described 
in [9]. The parameters for this model are PARJ(21) , PARJ(23) and PARJ(24) , which correspond, 
respectively, to the sigma of the main Gaussian, to the sigma and to the amplitude fractions for 
the second Gaussian. The result of all these improvements is the COMPASS tuning. 

In figure [1] the kinematic distributions and some hadronic distributions are shown, details can 
be found in the figures caption. In these figures, two MC simulations with different tunings are 
shown: the LEPTO default and the COMPASS tuning. In general the COMPASS tuning 
describes better the data, particularly for the hadronic variables pT and J2Pt which a 
remarkable agreement between MC and data is presented in the figure. The same applies for 
the hadron multiplicity (rightmost in figure [1]). In table [1] the values of the modified parameters 
are shown. 



Table 1. Modified JETSET MC parameters. 



Tuning 


PARJ(21) 


PAR J (23) 


PARJ(24) 


PARJ(41) 


PARJ(42) 


LEPTO Default 


0.36 


0.01 


2.0 


0.3 


0.58 


COMPASS 


0.34 


0.04 


2.8 


0.025 


0.075 




Figure 1. Data and MC comparison. Upper plots: (from left to right) distributions of pT of 
the leading hadron, pT of the sub-leading hadron, sum of the leading and sub-leading hadron 
p^, p of the leading hadron, p of the sub-leading hadron and hadron multiplicity. Lower plots: 
respective data over MC ratio. 



5. Systematic Studies 

The total systematic error is 6(AG/G) = 0.063. The contributions come from several sources 
which were studied in detail. For the neural network, the dependence on the internal structure 
was taken into account. For the MC several samples with different configurations (tuning, 
PDF, with and without PS, with and without Fl) were studied. An extensive study was 
performed searching for false asymmetries. The dependence of the asymmetry Af using different 
par ametrisat ions was studied. For equation (Ej) two approximations were used for x^: one with 
= const • xc, and another obtained from a second iteration in the neural network of xc 
given as input. Also the uncertainty for the input variables related to the polarisation states of 
the beam and target: P5, Pt and /, were taken into account. The major contributions come from 
the MC, S(AG/G)mc = 0.045, and the AG/G formula in equation (Ej), 5(AG/G)formuia = 0.035 . 

6. Results 

The gluon polarisation, AG/G, is calculated using equations ([2j) to ([5j). The result is AG/G = 
0.125 ± 0.060 ± 0.063 calculated at Xq^ = 0.091qq4 . In order to investigate xq dependence of 
AG/G the data is divided into three bins of the parametrised xq variable, i.e. divided into 
three independent samples. These results are given in table [2] and also presented in figure El 
In the same figure, other results from the COMPASS collaboration are also depicted, together 
with results from HERMES and SMC experiments. 



Table 2. Gluon polarisation results in bins of xq- 
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Figure 2. AG/G Results from COMPASS [3], SMC [E], HERMES [E] experiments. Also 
NLO QCD fits are shown [S]. 
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